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Abstract 

We report a measurement of the D^-D^ mixing parameter ycp using 23.4 fb~^ 
of data collected near the T(45') resonance with the Belle detector at KEKB. 
ycp is measured from the lifetime difference of mesons decaying into 
the K~TT~^ state and the CP even eigenstate K~K~^. We find ycp = 
(^—0.5 lb l.O^Qg^ X 10^^, where the first error is statistical and the second sys- 
tematic, corresponding to a 95% confidence interval —0.030 < ycp < 0.020. 

PACS numbers: 12.15.Ff,13.25.Ft,14.40.Lb 
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The search for D -D mixing provides an important window on physics beyond the Standard 
Model (SM). Since the predicted rate of mixing in the SM is very small Q, large mixing could 
indicate a contribution from non-SM processes. One measure of mixing effects is the lifetime 
difference between D^ decaying to the K~K~^ final state (a CP-even eigenstate) and the K~n~^ 
final state (which is not a CP eigenstate) 

_ t{K-tt+) 

= r{K-K+) - ^' 

where r = l/F and T is the effective decay rate obtained by fitting a single exponential to the 
measured decay distribution for each final state Q. We combine decays from D^ and D^, which we 
assume to be produced at equal rates in e'^e'^ collisions. The parameter ycp can be approximated 
as ucp ~ ycoscj) + xAsincf). Here x = (Mi — M2)/rav and y = (Fi — r2)/2rav, where Mi^2 and 
ri^2 are the masses and decay widths for the |-Di,2) = p\D^) ± q\D^) physical states of the neutral 
D meson system, and Tav = ^(Ti + r2). (j) is the phase of qA(D'^ K-K+)/pA{D^ K-K+), 
where A{D'^ (Z)°) K~K^) are the decay amplitudes, and A = (|pp - |gp)/(|p»P + \q\'^)- In the 
CP-conserving limit, ycp = U @- 

The FOCUS collaboration recently reported ycp = (3.42 ± 1.39 ± 0.74) x 10"^ Q. A common 
SM prediction is 0(10^^). Since non-SM processes are expected to enhance x but not y, 

such a large value of ycp would be difficult to interpret as a signal of new physics. Possible SM 
effects at the 10~^ level would be a more natural explanation |5|. We note however that limits on 
the mixing parameter x are weak (|x| < 0.03 ~ 0.06), and the parameter A is not constrained by 
existing measurements [||, so that a large ycp could also be accommodated if CP violating effects 
(A and/or (j)) were large. 

In this Letter, we present a new measurement of ycp with better statistical precision than 
the FOCUS result and largely independent systematic errors. The data sample for this analysis 
corresponds to an integrated luminosity of 20.9 fb~^ taken at the T(4S') resonance and of 2.5 fb~^ 
taken 60 MeV below the T(4S') resonance collected with the Belle detector at the KEKB col- 
lider. KEKB 1^ is an asymmetric energy electron-positron collider designed to produce boosted 
B mesons. The electron and positron beam energies are 8 GeV and 3.5 GeV, respectively: the 
resulting energy in the center-of-mass system (cms), 10.58 GeV, corresponds to the mass of the 
T(45) resonance. 

The Belle detector Q consists of a three-layer double-sided {rcf) and rz planes) silicon ver- 
tex detector (SVD), a 50-layer central drift chamber (CDC), an array of 1188 aerogel Cerenkov 
counters (ACC), 128 time-of-flight (TOF) scintillation counters, an electromagnetic calorimeter 
containing 8736 CsI(Tl) crystals and 14 layers of 4.7-cm-thick iron plates interleaved with a system 
of resistive plate counters (KLM). All subdetectors except the KLM are located inside a 3.4 m 
diameter superconducting solenoid which provides a 1.5 Tesla magnetic field. The impact param- 
eter resolutions for charged tracks are measured to be a'^y = (19)^ + (50/ (p/3 sin'^/^ O))^ //m^ in the 
plane perpendicular to the beam [z) axis and al = (36)^ + (42/(p/? sin^/^ d)f along the z axis, 
where /3 = pc/E, p and E are the momentum (GeV/c) and energy (GeV), and 9 is the polar angle 
from the z axis. The transverse momentum resolution is {cTpt/ptY = (0.0019pt)^ + (0.0030)^, 
where pT is in GeV/c. 

Candidate D^ mesons Q are reconstructed via Z?" K~t:~^ and D^ — > K~ decays. We 
require candidate charged tracks to be well reconstructed, and associated with at least two SVD 
hits in both the rcf) and rz planes. Charged pion and kaon candidates are selected based on 
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a particle identification likelihood value calculated using CDC energy loss measurements, flight 
times measured in the TOF, and the response of the ACC; for the chosen cuts, are identified 
with an efficiency of ~85% and a charged pion fake rate of ~10% for momenta up to 3.5 GeV/c. A 
candidate is required to have cms momentum greater than 2.5 GeV/c to eliminate secondary 
mesons originating from BB events. To suppress the background due to random combinations 
of two oppositely charged tracks, a cut is imposed on the decay angle 0d defined as the angle 
between the momentum vector of the candidate in the laboratory frame and that of the tt^ 
{K~^) in the rest frame for the K~tt^ {K~K~^) decay: cos^d > —0.85 for the K~tt^ decay 
and I cos^dI < 0.90 for the K~K^ decay, where the cos^d distribution is expected to be flat for 
signal. The two tracks forming the candidate are required to originate from a common vertex. 
In addition, the reconstructed flight path is required to be consistent with originating at the 
e'*"e~ interaction point (IP) proflle. Figure ^ shows invariant mass distributions for K~tt^ 
and — > K~K~^ candidates, superimposed with the result of a flt with two Gaussians (for signal) 
plus a linear function (for background). We find 214260 ± 562 K-it+ and 18306 ± 189 

K^K^ signal events, as determined from the area of the two Gaussians. The signal purity 
in the mass region within 3fT of the measured mean mass is 87% for K^-k^ (67% for K^), 
where a is the weighted average of the standard deviations of two Gaussians and is 6.5 MeV/c^ 
(5.4 MeV/c^ for K-K+). 

The decay vertex (a^dec) of the charm meson is determined using both tracks that form the charm 
meson candidate. The production vertex (ajpro) is obtained by extrapolating the fiight path to 
the IP. The center and size of the IP profile are determined for each KEKB fill |1C]. The size of the 
IP region is Ux = (80 — 120) /im, ay = {2 — 4) /um, and cr^ = (3 — 4) mm. The signed decay length 
in 3-dimensional space fill and the proper-time (t) are obtained from i = (a^dec — a^pro) • Pd/\Pd\ 
and t = lmD/\PD\^ respectively, where p£) and rriD are the momentum vector of the reconstructed 
charm meson and the world average value of the mass [12|. The selected mesons have a 



laboratory momentum of 3.6 GeV/c and a decay length of ~200 ^m on average. A Monte Carlo 
(MC) simulation study indicates that the decay and production vertex resolutions along the 
fiight direction are 110 /im (rms) and 70 //m (rms), respectively. We reject a small fraction (~ 3%) 
of the candidates by requiring that the uncertainty of the decay length measurement be less 
than 300 //m. 

In this analysis we extract the value of ycp by combining likelihood functions for K~tt^ 
and K^K^ decays, 

^ycp = ^K--K+ ■ Lx-K+i 
and expressing the lifetime for K~K^ as 

t{K-K+) = t{K-tt+)/{1 + ycp) 

in an unbinned maximum likelihood fit. This method allows us to properly estimate correlated 
systematic errors. The likelihood functions L = L^-^+ and Lj^-j^+ are given by |13] 

L{tsig,S, S'tail, /tail, •^BG, /tbg' '^'bG, 5't^il ' /taU ) 
roo 1 



n /siG / dt' e^SlG R{f - t'- S, 5tail, /tail) 

j L Jo rsiG 

rco 1 -t' 

+ (1-/sig) dt'if^^ e^BG 

Jo tbg 

+(1 - /rBG)'^(*')}-R(** - t'; crl, Sbg, Sf^i, /tan ) 
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with separate parameters for K~Tr^ and K~K^ , giving a total of 18 parameters to fit including 
ycp- The product is over the candidates. Here, /gjQ is the probability that the candidate is 
signal, calculated as a function of the invariant mass, and tsig is the signal lifetime. The function 
R represents the resolution of the proper time t*. The background proper-time distribution is 
modeled by a fraction /^^^ (~ 0.15 for K^tt^ and ~ 0.21 for K^K^) with the effective lifetime 
tbg (~ 391 fs for K^t:^ and ~ 497 fs for K^K^) and a fraction with zero lifetime represented by 
the Dirac delta function 6{t'). 

The resolution function i?, separately for the signal and the background, is parameterized as: 



where at is the proper-time error estimated candidate-by-candidate from the decay length error, 
and S and 5'taii are global scaling factors that modify at for the main and tail Gaussian distributions 
and /tail is the fraction of the tail part. The main component is due to the SVD vertex resolution, 
while the tail component is due to poorly reconstructed tracks {e.g. tracks affected by misassociation 
of SVD hits, wrong SVD hit clustering or large angle multiple-scattering). 

A simultaneous fit is performed to all the candidates contained in the mass region within 
40 MeV/c^ of the mean mass. This wide range includes both the signal (< ztScr) and the 
background-dominated (> itSo") regions. Figure ^ shows the results of the fit in the signal region. 
The solid lines show the fit and the dotted lines show the background contribution in the fit. 
Figure ^ shows the results of the fit in the background-dominated region and demonstrates that 
the background shape is well reproduced in the fit. The fit yields /tail ~ 0.17, S ~ 0.84, and 
5'taii ~ 1-75 resulting in an average proper-time resolution of 215 fs for signal, while ^ 0.06, 
5bg ~ 1-05 and Sf^^ ~ ^ fo'^ background. These values are found to be nearly identical for the 
K^TT^ and K~K^ decays. 

We obtain ycp = (-0.2±1.0) x 10"^ and t{K-t:+) = 416. 2± 1.1 fs from the fit. The fit results 
are corrected for small biases found in a large sample of MC events, where the reconstructed proper- 
time is smaller than the generated value. The difference diminishes when the effects of multiple 
scattering, decay in flight, or hadronic interaction are turned off in the MC simulation ||l^. The 
bias is found to be decay-mode dependent, —1.5 it 0.6 fs for K~tt'^ and —2.7 it 1.2 fs for 

K^K^, resulting in a correction for ycp of (—0.3 it 0.3) x 10^^. Here errors are due to the 
MC sample statistics and included as a systematic error in the final result. With this correction 
to the fit we obtain the result ycp = (—0.5 ± 1.0) x 10~^. 

The systematic uncertainties for ycp are listed in Table |[ Because ycp is measured as a ratio 
of two lifetimes, many correlated systematic uncertainties in the reconstructed decay length cancel; 
as a result, errors in decay vertex measurement, uncertainties in the IP profile and the uncertainty 
of the reconstructed momentum vector all make negligible contributions to the uncertainty of 
ycp- A MC simulation study shows that the fit yields a signal lifetime slightly longer than the 
generated one due to the presence of the background. This bias is common to the K~Tr^ and 
K~K~^ decays in sign and magnitude (~+3 fs) and, therefore, has no effect on the value of ycp- 

A difference in the signal purity between the K~tt^ and K^K^ decays may result in a bias in 
ycP- We estimate this uncertainty by repeating the analysis for samples obtained by varying 
requirements for particle identification and the fit quality of the decay vertex, both of which 



R{t] at, S, 5tail, /tail) = (1 - /tail) 



+ /tail 
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are very effective for suppressing background. The systematic uncertainty due to the background 
proper-time distribution is estimated by varying the mass range used in the fit from the nominal 
±40 MeV/c^ to ±35 MeV/c^ and to ±45 MeV/c^, and by comparing the results with diff'erent 
parametrizations. 

The systematic error originating from a correlation between the proper-time measurement 
and the measured mass is estimated by a MC simulation study. We estimate the systematic 
uncertainty due to candidates with large proper times {t > lOr^io) by varying the t range for 
the fit and taking the maximum excursion to be the systematic error. Also included are the effects 
of the systematic uncertainties due to the statistical uncertainty of the signal probability /sig and 
the error of the world average value of the mass |12|. The total systematic error is calculated 



by taking a quadratic sum of all contributions and is IqIoos- 

In summary, we have measured the D^-D^ mixing parameter ycp, using 23.4 fb^^ of data 
collected near the T(45) resonance, to be 



ycP = (-0.5 ± 1.0l°]) X 10- 



■2 



This corresponds to a 95% confidence interval —0.030 < ycp < 0.020. The result is consistent with 
zero and the Standard Model expectation that ycp is small. 
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Russian Federation; the National Science Council and the Ministry of Education of Taiwan; and 
the U.S. Department of Energy. 
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FIG. 1. The invariant mass distributions for (a) K^tt^ and (b) K~ candi- 

dates. The results of the fit with two Gaussians (signal) and a linear function (background) are 
superimposed. The dotted line indicates the background. 
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FIG. 2. Proper-time distributions and fit results for the decay modes (a) K~t:~^ and (b) 

— > K~K~^ in the mass signal region. The solid line is the result of the fit. The dotted line 
indicates the background contribution. 
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FIG. 3. Proper-time distributions and fit results for the decay modes (a) K vr^ and (b) 

K^K'^ in the mass background-dominated region. 



TABLE I. Summary of systematic errors on the ycp measurement. 



Source 


Systematic error 


(10-2) 


Reconstruction bias correction 




±0.3 


Decay vertex error 


ne^ 


eligible 


IP profile 


ne^ 


eligible 


momentum error 


ne^ 


eligible 


Particle identification 




±0.5 


Decay vertex quality 




+0.1 
-0.4 


Background t distribution 




+0.2 
-0.1 


mass - t correlation 




±0.3 


Large proper times 




±0.2 


Signal probability /^jq 




±0.1 


World average mass 




±0.1 


Total 




+0.7 

-0.8 
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